Non-inertial lateral migration of vesicles in bounded Poiseuille flow 
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Cross-streamline non-inertial migration of a vesicle in a bounded Poiseuille flow is investigated 
experimentally and numerically. The combined effects of the walls and of the curvature of the 
velocity profile induce a movement towards the center of the channel. A migration law (as a function 
of relevant structural and flow parameters) is proposed that is consistent with experimental and 
numerical results. This similarity law markedly differs from its analogue in unbounded geometry. 
The dependency on the reduced volume v and viscosity ratio A is also discussed. In particular, the 
migration velocity becomes non monotonous as a function of v beyond a certain A. 



Flow of confined soft entities, such as vesicles (closed 
quasi-inextensible lipid membranes), or blood cells in the 
circulatory system and in microfluidic devices, is a prob- 
lem of a paramount importance with both fundamen- 
tal and practical interests. While inertial effects can in- 
duce lateral migration of any flowing body in a chan- 
nel [H, 0, 3] , the ability of these soft entities to adapt their 
shapes under non-equilibrium conditions gives them the 
possibility to migrate transversally even at low Reynolds 
number. Transverse migrations induce non uniform lat- 
eral distributions of the suspended entities, which has 
important consequences on the rheology of a confined 
suspension (e.g. the Fahraeus-Lindquist effect in blood 
vessels [H), or should impact on transport efficiency in 
the various sorting microfluidic devices that are now be- 
ing developed @. 

Despite the considerable interest for transverse migra- 
tion in many circumstances, there is, to our knowledge, 
yet no quantitative law that would allow one to relate the 
lateral migration velocity of a deformable entity flowing 
in a channel, even isolated, with its position, its mechan- 
ical properties, and the flow parameters. We propose a 
law for the case of a single vesicle placed in a bounded 
Poiseuille flow, from experiments in microfluidic devices 
as well as simulations based on the boundary integral 
method. 

The behavior of vesicles under unbounded shear flow 
has been the subject of several theoretical @, 0, H, H El 
and experimental [ll|, E2, El studies. When the vis- 
cosity ratio between the inner and the outer fluids is 
small, vesicles perform a tank treading dynamics where 
the orientation of the main axis of the vesicle is con- 
stant and the membrane undergoes a tank treading mo- 
tion. In a bounded Poiseuille flow, tank-treading vesi- 
cles experience a transverse force and reach the center 
where they assume a steady shape. The latter stage 
has been described in several papers [3, EH, El- Mi- 
gration in Poiseuille flows has been also reported on 
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capsules EzLGJ, El, red blood cells [H, El, and 
drops |20l. l2l|. The latter seem to have a very differ- 
ent behavior: depending on the viscosity ratio and the 
confinement, the reported equilibrium positions are not 
always on the centerline. In addition, a drop interface 
and a vesicle membrane are mechanically different: a 
vesicle has a constant surface area and its equilibrium 
shape is not a sphere in general. As a consequence of 
these intrinsic differences the vesicle shape under flow is 
described (in the co- moving frame), to leading order, by 
a nonlinear shape equation while a drop obeys a linear 
equation Q. 

Lateral migration has two distinct sources, (i) A wall- 
induced lift force [12, [H 0, In agreement with 
the numerical and theoretical studies for simple shear 
flows [H, [24[ , we showed recently that the migration ve- 
locity decreases like 1/y 2 , where y is the distance to the 
wall of the center of mass of the vesicle [11 . (ii) The non- 
constant shear rate in a parabolic velocity profile (even 
unbounded) leads to a subtle interplay between the gra- 
dient of shear and the shape (27[, resulting in migration 
towards the center with a constant drift velocity except 
near the center-line. In a realistic channel, both effects 
coexist and we shall see that this leads to a new and 
nontrivial non-inertial migration law. 

The considered microfluidic channel is straight and has 
a rectangular cross section. The flow direction is Ox, and 
we investigate lateral migration along Oy. Let 2uu denote 
the channel width in the y-direction, and vq the imposed 
flow velocity at the center of the channel in the absence of 
vesicle. The two walls are located at y = and y = 2w. 

A vesicle is characterized by two geometrical parame- 
ters: its effective radius Rq, determined from its constant 
volume V by J?o = (3V/47T) 1 / 3 , and its reduced volume 
v = V/[47r(5/47r) 3 / 2 /3] {S is the constant area of the 
vesicle) characterizing vesicle deflation. Volumes are cal- 
culated at lift-off (see below) by assuming axisymmctric 
shape about the vesicle's main axis. The viscosity ra- 
tio is defined as A = rji n /rj out , where r\i n ,r\ out denote 
the inner and the outer viscosities. Relevant space and 
time scales are the vesicle radius Ro and the characteris- 
tic time needed by the vesicle to relax to its equilibrium 
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FIG. 1: (color online) (a): Experimental time evolution of 
the lateral position y for 3 vesicles of similar size (19±1 /xm) 
and imposed flow velocity vo — 920/im.s -1 but different v 
and A. The dashed line indicates the center- line (which will 
be reached much later by two of the vesicles). The solid line 
shows the y(t) curve obtained from the fit by the numerical 
solution of Eq. [T] (b) : Evolution of y — yo versus t for vesicles 
with A = 1.1 and different v. For each v interval, the different 
symbols correspond to different vesicles, whose w and vq vary 
in the explored intervals, that are 1.8 < w < 9.6 and 400 < 
vo < 14800. The corresponding curves clearly collapse on a 
single one. For clarity, the sets of curves corresponding to 
different v intervals are switched vertically with an increment 
of 0.5. 



shape (in the absence of imposed flow), which is given 
by t = r] ou tR%l n, where n ~ 20fcsT is the membrane's 
bending rigidity (typically, r ~ 10 s). 

Therefore dynamics depends a priori on the four di- 
mensionless parameters (w = w/Rq, v = vqt/Rq, u, A). 
We consider first the case A ~ 1. The strategy is to vary 
w and the imposed velocity vq and investigate the migra- 
tion law y(i) — y(t/r)/Ro for each value of v. Then we 
discuss the influence of v and A. 

In the experiments, we used straight channels of height 
ho = 66.6 fj,m (in the direction of gravity z) and width 
2w (rectangular cross section) between 70 and 140 /im. 
The walls of the channels are made of PDMS glued to 
a glass slide. The flow is induced by gravity, by con- 
necting the inlet and the outlet to reservoirs at different 
heights. Vesicles are prepared following the electroforma- 
tion method. They are made of a dioleoylphosphatidyl- 
cholinc (DOPC) lipid bilayer enclosing an inner solution 
of sugar (sucrose or glucose) in water or in a 1:4 glycerol- 
water (w:w) mixture. Samples are diluted in a slightly 
hyperosmotic outer solution of the same type, in order to 
deflate them by osmosis. Dextran can be added to one of 
the solutions to modify the viscosity ratio A. Vesicle size 



i?o lies in the range 7-37 [im while vq varies between 200 
and 1100 /jm.s -1 . Note that for our solutions of viscos- 
ity and density close to the one of water, the Reynolds 
number R e = pvoRo/Vout is always lower than 4.10 -2 . 

A particular design of the upstream channel creates 
an initial condition where incoming vesicles touch the 
y = wall in the observation area and start to be lifted 
away from it. In particular, they have already devel- 
oped a nearl y el lipsoidal shape tilted with respect to the 
wall [H, [24], HI] . The 2D fluid velocity profile in the xy 
plane where their center of mass lies is nearly parabolic, 
since the rectangular cross section of the 3D channel 
obeys 2w/ho < 3 (28|. Moreover, we wait for the flow to 
be established for a long time, resulting in preliminary 
centering of the vesicles in the z direction. The imposed 
velocity profile is thus written as v"^(r) = c (yw — y 2 /2), 
where c = 2vq/w 2 is the curvature. A vesicle is tracked 
along its trajectory with a phase contrast microscope, 
and the position y of its center of mass is determined by 
image processing, starting with y(t = 0) = yo, where yo 
is the position just before lift-off, which is close to Ro- 

The evolution with time of the y position of two vesi- 
cles with A = 1.1 is shown in Fig. QJa). The vesicles 
quickly move away from the wall, then the migration 
velocity decreases to zero as they approach the center- 
line. Along their trajectory, they continuously deform 
from a tilted ellipsoid to a symmetric bullet-like shape. 
For a given v, the function y(t) depends a priori on 
both parameters w and «o- In order to determine this 
functional dependence, which is not known a priori, we 
rescale the time variable. The choice of a relevant time 
scale is not obvious, however. Indeed, while the inverse 
of the shear rate yields a natural scale, this is not an ad- 
equate choice since the shear rate is not constant along 
the trajectory. The trick is to rescale each infinitesimal 
time step dt around the time t by the local shear rate 
j(y) = dv°° / dy — c{w — y) of the unperturbed flow at 
the position y(t). The new dimensionless time-like pa- 
rameter is then obtained by integrating the rescaled time 
steps: t = J Q t j{y)dt' = c f Q [w — y(t')\dt'. t accounts for 
the history of the shear rates experienced by the vesicle 
along its trajectory. The raw data (not shown) for the 
migration velocity spreads over more than a decade in the 
parameters space. Interestingly, as shown in Fig. QJb), 
all experimental curves y(t) for a given v (or rather, a 
tiny interval around it) collapse, whatever the values of 
w and vo within the explored range. A log-log plot (not 
shown) of each master curve y(t)—yo is linear, a clear sig- 
nature of a power law behavior y(t) — yo = (3t a , where the 
dimensionless parameters a and are thus independent 
from w and vq. The lateral migration velocity v m = y 
as a function of the position y and (w,Vq) is then easily 
extracted: 
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that is: y = £ 



Ro +L i(y) 
(y - yo) 5 ' 



(i) 



This law constitutes the central result of our find- 
ing. In the range 0.970 < v < 0.975 we find for 
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FIG. 2: (color online) Simulations, (a): migration velocity 
v m versus lateral position for different reduced volumes. Solid 
lines show the fits to Eq.[T] (b): log- log plot of v m versus y—1 
for different w and vo and v = 0.97. 

instance £ = a/3 1 /" = 1.2 x 1CT 2 ± 0.2 x lfr 2 and 
S = 1/a — l = 1±0.1. The error bars for these coefficients 
are mainly due to the uncertainties on the measure of u, 
of order ±0.005, since, as we shall see, the velocity de- 
pends on the reduced volume. Note that the differential 
equation [T] has no analytical solution but can be easily 
solved numerically and the result used to fit the raw data 
y(t) without rescaling procedure (Fig. QJa)). 

In the simulations, we studied 2D neutrally buoyant 
vesicles (the 2D geometry captures the essential features) 
having no viscosity contrast (A = 1). The fluids flows 
inside and outside the vesicle are governed by the Stokes 
equations (r m £ dfl is a membrane point): 

-Vp(r)+??V 2 v(r) = -5(r-r m )f(r), V-v(r)=0, (2) 

where p is the pressure, v is the velocity and f the mem- 
brane force, given by Eq. 8 of Ref. 27]. Thanks to the lin- 
earity of Eqs. [2J we solve them using a boundary integral 
method adapted to vesicle problems @, 0, H3, H(| • 
The membrane velocity is then given by the following 
integral equation, that we solve numerically: 

Vi(r m ) — —— (j) G^(r m ,r')^(r')rf S (r')+«r(r m ), (3) 
47r7 7 Jan 

where Gjj is the Green's function for a fluid bounded by 
a steady infinite plane wall located at y = 0: 

G^(r,r') = Gyfr-rO-G^r-rJ) 

+ 2rfGg(r - r^) - 2r' y Gff(r r' r ). (4) 

Gjj(r) = — Sijhxr + ^J 2 - is the Green's function for an 
unbounded fluid, or Stokeslet, r'j — (r' x , —r' y ) is the image 
of r' with respect to the wall. The function 

G§(r) = (^-2^) (5) 

is the Stokeslet doublet, and the source doublet is 

Gf(r) = r,Gg(r) + (6 jx - s j^^ . (6) 
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FIG. 3: (color online) Reduced migration velocity v m of a 
vesicle at position y — 4 versus its reduced volume. For read- 
ability, the simulations data are uniformly rescaled by a factor 
~ 0.1. 



The evolution of the vesicle's shape and location are ob- 
tained by updating every membrane point using a Eulcr 
scheme: r m (i + At) = v(r m , i)df + r m (i). 

Note that we only consider one wall at y = 0. Provided 
w > 8, the vesicle reaches and stays on the y = w line at 
long times, even without the symmetric wall at y = 2w. 
This seems to indicate that for w > 8, migration forces 
due to the curvature dominate over wall effects near the 
center-line. 

The variations of v m with y are shown in Fig. HJa) for 
four different reduced volumes and given width and flow 
velocity. They are well described by Eq. [1] The ade- 
quation to this law is confirmed by considering varying 
8 < w < 12 and 600 < v < 800. Following Eq. Q] it 
is convenient to rescale the migration velocity in such a 
way that it should not depend either on w, nor on #o : 
v m = v m /j(y) ih) 5 - The variations of v m with 

y — 1 are shown in Fig. [IJb) for v = 0.97. The collapse 
on a single line (in log-log scale) is in accordance with 
Eq.[Q We find 5 ~ 0.8 ±0.1 and f ~ 0.1. The agreement 
between experiments and simulations regarding the expo- 
nent S is quite satisfactory. However, numerical studies 
overestimate the amplitude £. This is attributed to the 
2D character (actually a translationally invariant form 
in the z direction), causing an enhancement of the lift 
force. While 2D simulations have captured some inter- 
esting facts, a 3D simulation is necessary before drawing 
general conclusive answers, and we plan to investigate 
this problem in the future. 

We now discuss the dependencies on v and A; for all 
the values of v and A explored here, the experimental 
and numerical curves are still very well fitted by the law 
given by the resolution of Eq. [T] (see Fig. HJa)). Values 
for v m at y = 4 are reported on Fig. [3] For A ~ 1, in 
the ranges 0.83 < v < 1 and 0.89 < v < 0.99 offered 
respectively by the experiments and the simulations, we 
find an increasing migration velocity with decreasing v 
(as seen on Fig.[2ja)). 

Results for higher A are available from the experiments. 
The new feature is the non monotonous behavior of v m 
as a function of v. When v is decreased from 1, the 
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migration velocity first increases (as for A = 1.1), then 
reaches a maximum and decreases back to a very low 
value (Fig. El). 

This non-monotonous behavior of the velocity can be 
understood on a general physical basis. A spherical vesi- 
cle (y — 1), does not migrate owing to the fore-aft sym- 
metry. On the other hand, when v decreases, the vesicle 
finally switches from a tank-treading to a tumbling mo- 
tion [la| . In the latter regime, no global migration should 
occur either because the averaging over the different ori- 
entations of the vesicle during one rotation period leads 
to an almost symmetrical configuration. From these con- 
siderations we infer a maximal velocity at a given value 
of v. Note that when A ~ 1, no tumbling motion occurs 
whatever v, so that a monotonous evolution of the veloc- 
ity with v is observed. The evolution with v and A is in 
qualitative agreement with the predictions made by Olla 
for the migration velocity of a (shape preserving ellip- 
soidal) vesicle in the case of a simple shear flow bounded 
by a wall, although the scaling is different (22J. 

In conclusion, our experiments and simulations yield 
a similarity law for the lateral migration velocity of a 
vesicle in a bounded Poiseuille flow as a function of its 
distance to the walls and to the center-line, its effective 



radius, the channel's width and the flow velocity. We 
showed that the effects of the walls and of the curvature 
of the velocity field are coupled in a non linear manner: 
curvature not only induces migration [27| but also affects 
the shape and orientation, which affects the lift force. 
The law v m ~ j{y)/y markedly differs from what the 
naive extrapolation of the results for a vesicle near a wall 
and in a linear shear flow would give: v m ~ A f{y)/y 2 - 

Deflating a spherical vesicle increases its deformability, 
thus its asymmetry under shear, and leads to higher 
migration velocities. However, beyond a given viscosity 
ratio, the tank-treading to tumbling transition is ap- 
proached when the deflation increases, and the migration 
velocity undergoes a decline which can be understood 
on the ground of general symmetry considerations. 
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